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Abstract 


We continue our study of entanglement entropy in the holographic super- 
conducting phase transitions. In this paper we consider the holographic p-wave 
superconductor /insulator model, where as the back reaction increases, the tran- 
sition is changed from second order to first order. We find that unlike the s-wave 
case, there is no additional first order transition in the superconducting phase. 
We calculate the entanglement entropy for two strip geometries. One is par- 
allel to the super current, and the other is orthogonal to the super current. 
In both cases, we find that the entanglement entropy monotonically increases 
with respect to the chemical potential. 
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1 Introduction 


The entanglement entropy nowadays plays an important role in understanding some charac- 
terizations in quantum field theory and many-body physics (see, for example, Refs. [1, 2, 3]). 
The entanglement entropy, as the name suggests, measures how a given quantum system is 
entangled or strongly correlated. Further, it can also be used to distinguish different phases 
and corresponding phase transitions and is considered as a useful tool for keeping track of 
the degrees of freedom of strongly coupled systems. For a given system, the entanglement 
entropy of one subsystem with its complement is defined as the von Neumann entropy. 
However, the calculation of entanglement entropy is found to be very difficult except for 
the case in 1 + 1 dimensions. 

In the framework of AdS/CFT correspondence |4, 5, 6], a holographic method to cal- 
culate the entanglement entropy has been proposed in Ref. [7]. Following Ref. [7], for a 
conformal field theory (CFT) which has a dual gravitational configuration living in one 
higher dimension, the entanglement entropy of the CFT in a subsystem A with its com- 
plement is given by finding the minimal area surface y4 extended into the bulk with the 
same boundary OA of A (see Refs. [8, 9] for reviews). Namely, the entanglement entropy 
of A with its complement is given by the “area law” 


QT 
SA = 2 Area(ya), (1) 


where k? = 87G is related to the gravitational constant in the bulk. 

Holographic superconductor model was first constructed in Refs. [10, 11] where the 
model is a s-wave one since the condensed field is a scalar field dual to a scalar operator in 
the field theory side. The holographic approach was also generalized to the p-wave case [12] 
and d-wave case [13, 14]. It is interesting to note that the back reaction of matter field on 
the background geometry will change the order of the holographic p-wave superconductor 
phase transition from second order to first order [15]. On the other hand, the effect of higher 
derivatives on the holographic p-wave model was studied in Refs. [16, 17, 18]. And near crit- 
ical point, an analytic study on the holographic p-wave superconductor/conductor phase 


transition has been made in the probe limit [19, 20]. In Ref. [21], a holographic s-wave su- 
perconductor/insulator phase transition model was built at zero temperature. Combining 
the conductor/superconductor phase transition with the insulator/superconductor phase 
transition, Horowitz and Way [22] described a complete phase diagram for a holographic 
s-wave superconductor/conductor/insulator system. In addition, the holographic p-wave 
superconductor/insulator phase transition was discussed in Ref. [23], and Ref. [24] carried 
out an analytic study for the holographic superconductor /insulator phase transition. While 
various properties of holographic superconductor models have been intensively investigated 
in the literature, in this paper we are interested in the behavior of entanglement entropy in 
the holographic superconductor models. Refs. [25, 26, 27] studied the behavior of entangle- 
ment entropy in holographic conductor/superconductor phase transition models, including 
s-wave and p-wave cases. It showed that the entanglement entropy decreases as one lowers 
temperature, indicating the degrees of freedom are reduced at lower temperature, and that 
the entanglement entropy can tell us not only the occurrence of the phase transition, but 
also the order of the phase transition. It indicates that entanglement entropy is indeed 
a good probe to phase transition. In Refs. [28, 29] we investigated the behavior of en- 
tanglement entropy in the holographic s-wave superconductor/insulator model, and found 
that the entanglement entropy as a function of chemical potential is not monotonic in the 
superconducting phase: at the beginning of the transition, the entropy first increases and 
arrives at its maximum at some chemical potential, and then decreases monotonically. 

In this paper we continue our study of the entanglement entropy in holographic su- 
perconductor phase transitions. We here focus on the holographic p-wave superconduc- 
tor/insulator model, paying main attention on whether those properties of entanglement 
entropy found in the previous studies are universal or not, in particular, on the non- 
monotonic behavior of the entropy in the superconductor/insulator phase transition. This 
will be helpful to further understand the superconductor/insulator transition. Further- 
more, there are also other motivations. It is well known that there is no thermal phase 
transition in (2+1)-dimensional field theories due to the Mermin-Wagner theorem or Cole- 
man theorem, which states that in quantum field theory, continuous symmetries cannot be 
spontaneously broken at finite temperature in systems with sufficiently short-range inter- 
actions in spatial dimensions d < 2, which at least tells us that the phase transition can 
be affected by the space-time dimension of the system. However, one indeed observes the 
superconducting phase transition in the gravity dual of (3 + 1)-dimensional AdS black hole 
backgrounds [10]. This might be caused by the suppression of the large fluctuations in the 
large N limit, which is supposed to be one of conditions for the validness of the AdS/CFT 
correspondence. In the framework of AdS/CFT correspondence, we can study the strongly 
coupled systems simply through their dual weakly interacting gravity theories. So a natural 
question arises that how the order of the phase transition of the strongly coupled systems 
can be affected by the dimension of space-time. Indeed, it is observed in Ref. [27] that in 
contrast to the (3 + 1)-dimensional holographic p-wave superconductor/conductor phase 
transition, the phase transition in the (2 + 1)-dimensional case is always second order, 
regardless of the strength of the back reaction. To further stress the issue and enrich our 
understanding of the holographic superconductor phase transitions, therefore it would be 


helpful to study the full back reacted holographic p-wave superconductor/insulator model. 

The p-wave superconducting phase is characterized by the condensation of the vec- 
tor “hair” [30, 12], while the normal insulator phase is described by a pure AdS soliton 
solution. As one increases the chemical potential, the pure AdS soliton background will 
become unstable to develop the vector “hair”. The emergence of the “hair” induces the 
symmetry breaking and gives a finite vacuum expectation value of the dual current oper- 
ator in the field theory side, which plays the role of order parameter in the holographic 
phase transition. We will choose the condensate of the current operator along one spatial 
direction. Thus, the U(1) symmetry as well as the spatial rotational symmetry are broken 
in the superconducting phase. The order of the transition from the insulator phase to 
superconducting phase can be either second order or first order, depending on the strength 
of the back reaction of matter field on the bulk geometry. For the holographic s-wave 
superconductor /insulator case, one can find a new first order phase transition inside the 
superconducting phase [22, 29]. We will see that this does not happen in the p-wave model. 
In view of the anisotropy in the superconducting phase, we will consider two strip geome- 
tries with a finite width along directions parallel to the super current and orthogonal to 
the super current, respectively. The value of entanglement entropy calculated in former 
setup is larger than the one in the latter case for given parameters. In addition, no matter 
of the order of the transition, we find the behavior of the entanglement entropy as a func- 
tion of chemical potential is monotonic. More precisely, the entanglement entropy always 
increases as the increase of chemical potential. Such result is quite different from the one 
in the s-wave case. 

This paper is organized as follows. In section 2, we introduce the holographic model 
and study the thermodynamics of the transition. Since the system is anisotropic in the 
superconducting phase, we calculate, in section 3, the entanglement entropy for strip ge- 
ometry along spatial direction x and y, respectively. The conclusion and some discussions 
are included in section 4. Some numerical details are given in appendix A. 


2 Holographic P-wave superconductor /insulator model 


2.1 Gravity background 


Let’s begin with the (4+ 1)-dimensional SU (2) Einstein- Yang-Mills theory with a negative 
cosmological constant [12] 
Fe) - PF (2) 


1 
5 
s= fa sv -glz a(R + Er 


where g is the Yang-Mills coupling constant and L is the AdS radius. The SU(2) Yang-Mills 
field strength is given by t 


12 1 


Fe, = ð, A; — 8, At, + eM A? AS, (3) 


‘u,v = (t,r, x,y,z) denote the indices of space-time and a,b,c = (1,2,3) are the indices of the SU(2) 
group generators T° = o° /2i where o° are Pauli matrices. 
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where e% is the totally antisymmetric tensor with e! = +1. The gauge field is given by 


A = Alir*dx". Here we define a parameter a = «/g which measures the strength of the 
back reaction of the Yang-Mills field on the background geometry. 
Our ansatz for the metric and Yang-Mills fields are given by ? 


ds? = — —_ + 7?(—f(r)dt? + h(r)dx? + dy? + g(r)e dr’), (4) 


A= o(r)r3dt + w(r)r'dz. (5) 


g(r) vanishes at the tip r = ro of the soliton. Further, in order to obtain a smooth geometry 
at the tip ro, 7 should be made with an identification 


An Lex 
n~t, T= 
r69'(To) 


This gives a dual picture of the boundary theory with a mass gap, which is reminiscent of 
an insulator phase [21]. 

The independent equations of motion in terms of the above ansatz can be written as 
follows: 


(6) 
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X-F Gg h m rw fh? BS fgh 
where the prime denotes the derivative with respect to r and the tip r = rọ is determined 
by the condition g(ro) = 0. 

In order to match the asymptotical AdS boundary, the matter and metric fields near 
the AdS boundary r —> oo should behave as 


Q 
b=by- g ts w = wo + -z +t , 
8 
fa g4 h4 X4 (8) 
ae ge 5 TRI T ; h=1+-7t 3 =n j 


? Our ansatz adopted here is very similar with the one in Ref.[23]. In the case a — 0, the probe limit 
can work very well. Therefore one can recover the results obtained in the probe limit [23, 24] when one 
takes a = 0. 


where wọ is regarded as the source of the dual current operator. To break the U(1) gauge 
symmetry spontaneously, we impose wọ = 0. In addition, we choose the condensate along 
the x direction. Thus, the rotational symmetry in x — y plane is also destroyed in su- 
perconducting phase. According to the AdS/CFT dictionary, we can obtain the current 
condensate (Jt) = 20° Wa, chemical potential u = ¢@ 9 and charge density p = 20° by by 
reading off the coefficients w2, œo and @2 in (8), respectively. Furthermore, here f4, g4, ha 
and y4 are all constants. Let us note that there exists an analytic solution for the equations 
of motion (7), which is just the so-called AdS soliton solution 


o(r)=H, w(r)=x(r)=0, f(r)=Alr)=1, g(r) =1- J. (9) 


Here u is an arbitrary constant, but less than its critical value, which will be showed shortly, 
and r = ro stands for the tip of the soliton. The AdS soliton solution has vanishing w, 
thus corresponds to the normal insulator phase. In order to compare physical quantities 
obtained from different solutions, the boundary geometry must be the same, in particular, 
the 7 coordinate must have the same compactification length I. Therefore, we scale all of 
I’ for each solution to be mL. We work in unites where L = 1 hereafter, so when L appears 
in later sections it is set to be one. Numerical details for solving the equations of motion 
(7) can be found in appendix (A). 
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Figure 1: The condensate of current operator (J%) (left plot) versus chemical potential u 
for a = 0.2. T is scaled to be tL. The critical chemical potential in this case is ue ~ 2.265. 
It is a typical second order phase transition which can be seen clearly from right plot, 
which presents the grand potential of the soliton with vector “hair” (solid purple) and of 
the soliton without vector “hair” (dashed black), with respect to chemical potential. 


2.2 Thermodynamics and phase transition 


To find the thermodynamically favored phases and to determine the critical point of phase 
transition, we should calculate the grand potential function Q of the system since we are 
working in grand canonical ensemble in this paper. In gauge/gravity duality the grand 
potential Q of the boundary thermal state is identified with temperature times the on-shell 


bulk action in Euclidean signature. The Euclidean action must include the Gibbons- 
Hawking boundary term for a well-defined Dirichlet variational principle and further a 
surface counterterm for removing divergence 


1 ee E. 6 
S Euclidean = - | boa gaR+ 7) _ age r J+ aj davh(—2K +7), (10) 


where h is the induced metric on the boundary, and KX is the trace of the extrinsic curvature. 
By using of the equations of motion (7) and the asymptotical expansion of matter and 
metric functions near the boundary, the grand potential Q turns out to be ? 


21K 
V3 = Ya, (11) 


where V3 = f dxdydn. Since we have scaled I to be mL, g4 = —1 for the pure AdS soliton 
solution, namely, in the normal insulator phase. 

The condensation of vector operator Je and the grand potential as a function of chemical 
potential for a = 0.2 and a = 0.6, as two typical examples, are presented in Figure (1) and 
Figure (2), respectively. Above the critical chemical potential ue, the solution with non- 
vanishing vector “hair” is thermodynamically favored over the pure AdS soliton solution. 
The emergency of the vector “hair” breaks the U(1) gauge symmetry as well as rotational 
symmetry and represents the superconducting phase. In contrast to the case in Figure (1) 
where (J*) rises continuously from zero at fic, the case in Figure (2) has a jump in the 
(J2) at uc indicating the first order of the transition. The corresponding charge density is 
drawn in Figure (3), which behaves qualitatively similar as the condensation. The order 
of the phase transition can be seen more clearly from the grand potential Q with respect 
to chemical potential. For the first case, the grand potential decreases monotonically as 
one increases the chemical potential, indicating a continuous phase transition. For the 
second case, the grand potential with respect to chemical potential develops a “swallow 
tail”, which is a typical feature in first order transition. Thus, as expected, we find that 
the order of the phase transition depends on the parameter a, i.e., the strength of the back 
reaction. The critical value is a, ~ 0.538 + 0.002, beyond which the order of the phase 
transition changes from second order to first order. However, unlike the s-wave case, there 
is no additional first order transition occurring in the superconducting phase [22, 29]. 

In addition, it is worth pointing out that except for the pure AdS soliton solution (in- 
sulator phase) and hairy soliton solution (superconducting phase) mentioned above, the 
action (2) admits another two solutions: the black hole solution without hair (conductor 
phase) and black hole solution with vector hair (another superconducting phase). A rough 
phase diagram has been drawn in Ref.|23]. The details depend on the strength of the 


3The soliton background has no horizon and the associated Hawking temperature vanishes. But one 
can introduce an arbitrary inverse temperature 1/T as the period of the Euclidean time coordinate. Since 
the soliton solution is static, the integration over the Euclidean time in the Euclidean action just gives 
the factor 1/T, which cancels the temperature factor in the grand potential and leads to a finite grand 
potential. 
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Figure 2: The condensate of current operator (J) (left plot) versus chemical potential y 
for a = 0.6. T is scaled to be tL. The right plot shows the grand potential of the soliton 
with vector “hair” (solid red) and the soliton without vector “hair” (dashed black) with 
respect to chemical potential. It is a typical first order phase transition and the critical 
chemical potential denoted by vertical dotted line is ue ~ 2.157. 


back reaction. The dependence of the order of phase transition on the back reaction is 
mentioned there. In this paper we just focus on the case with zero temperature insul- 
tor/superconductor phase transition. In this case, the black hole phase does not appear. 
In particular, in the zero temperature case, for any values of chemical potential beyond 
the critical one, the superconducting phase always exists. 


3 Entanglement entropy 


Let us begin to study the behavior of entanglement entropy in this p-wave model. Since 
the system is anisotropic in the superconducting phase, we will consider a straight belt 
geometry with a finite width Z along spatial direction x and y, respectively. 


3.1 Strip with finite width along x direction 


We first consider a straight belt geometry with a finite width £ along the x direction. 
More specifically, the subsystem A sits on the slice r = + where e€ > 0 is the UV cutoff, 
and extends in y and ņ directions. The holographic dual surface y4 is defined as a three- 
dimensional surface 


i=0, r=r(z) -2 <y < Ž (R> œ), 0<n<Tr. (12) 


R is the regularized length in y direction. For this boundary subsystem, there exist two 
kinds of dual surfaces in the bulk. One is called connected configuration and the other is 
disconnected configuration. We first focus on the connected configuration which is smooth 
at the turning point r = rą: The holographic surface y4 starts from x = £ at r = E, 


extends into the bulk until it reaches r = r,, then returns back to the AdS boundary r = 4 
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Figure 3: The charge density as a function of chemical potential for a = 0.2 with ue ~ 2.265 
(left plot) and a = 0.6 (right plot), respectively. T is scaled to be mL. The critical chemical 
potential denoted by vertical dotted line is ue ~ 2.157 for the right plot. 


at x = -£ For such an embedding, following the formula (1), the entanglement entropy 
of the subsystem is given by 


4rL E h(r)e-X) Qn L 1 
eet = — T RT r fate) (r Je a NEN ee = T RI(5 ae Gen), 
E 


(13) 
where the UV divergent part 1/e? has been separated from the total entropy and S°” is a 
finite part. The width £ of the subsystem A and r, are connected by the relation 


= = L 
5f a= | eT (14) 
2 rs r re 72y/g(r)h(r)(—< glr)hlrje- x0 1) 


r$g(ra)h(rs)e™ x0 


— r°g(r.)h(re)e7X*) 


On the other hand, the disconnected configuration consists of two separated surfaces that 
are located at £z = +£ and both start from the AdS boundary r = 1/e, extend into the 
bulk until they reach the tip of the soliton solution r = rọ. The entanglement entropy for 
this disconnected configuration is independent of /, and given by 


, An L z ak 27L 
Sy alta = is rr | re” mr — T TERNG i gerseen; (15) 
K i 
where Ss” is the finite part. 


Once have numerically solved the metric functions, we can now calculate the entan- 
glement entropy. Following the above discussion we must calculate the entropies of the 
connected configuration (13) and the disconnected one (15), respectively. One must take 
care of the UV divergence for the calculations of each entropy, but the difference of en- 
tropies is insensitive to it. Therefore, what we are interested in is the finite part of the 
entanglement entropy, i.e., S@*°°" and S°” in our discussion. We find that the entan- 
glement entropy with respect to strip width ¢ behaves quite similar for different choice 
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Figure 4: The typical configuration for the entanglement entropy as a function of strip 
width /. The entanglement entropy plotted in this figure and hereafter is the finite part 
with the UV divergent part subtracted off. The dashed purple and solid blue curves come 
from the connected configuration, while the solid red one comes from the disconnected 
configuration. The lowest curve is physically favored compared with others. 


of parameters, i.e., œ and u, which is schematically drawn in Figure (4). We stress here 
that the entanglement entropy drawn in Figure (4) and hereafter is the finite part with 
the UV divergent part subtracted off. + We find that there exists a maximal width lmaz- 
For l > max, the connected configuration does not exist and the physical solution comes 
from the trivial disconnected configuration where the entropy is independent of @. On the 
contrary, there are three different branches for the entanglement entropy when £ < bmaz- 
Two of them come from connected configurations and the third one corresponds to the dis- 
connected one. The physical entropy can be found by always choosing the lowest branch 
since it has minimal value. As we can see from Figure (4), the lower connected branch 
is favored for sufficiently small 4. Hence, there is a critical value Zeriticaı below which the 
connected branch is physically favored. Thus, as we change £, a phase transition must 
occur at Ceritical, Which is just the so called “confinement /deconfinement” phase transition 
discussed intensively, for example, in Refs. [31, 32, 33]. To be precise, for 0 < Ceritical; 
the entanglement entropy comes from the connected configuration and exhibits non-trivial 
dependence on ¢, which describes a “deconfining” phase. For € > oritical, the entropy 
is dominated by the disconnected configuration and is ¢ independent, which indicates a 
“confining” phase. 

Thus we have totally four phases in the boundary field theory side, i.e., the confin- 
ing/deconfining phases for the insulator and superconductor [28], respectively, which are 


4We omit the superscript used to distinguish the connected embedding and disconnected one for con- 
venience. This will not make any confusion since it is easy to distinguish these two cases from the context. 
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Figure 5: The entanglement entropy as a function of chemical potential for strip with finite 
length along x direction. We choose a = 0.2 at fixed £ = 0.5L (within “deconfining phase”, 
left plot) and 2 — oo (within “confining phase”, right plot), respectively. The transition 
from the insulator to superconductor is second order here. In both plots, the solid curves 
come from superconducting phase, while the dashed black lines come from insulator phase. 
The critical chemical potential here is ue ~ 2.265. The physical curve is determined by 
choosing the dashed line below ue and solid curve above ue. 


characterized by the chemical potential u and strip width @. In particular, the strip width 
controls the “confinement /deconfinement” phase transition. ° 

It is instructive to study how the entropy of the subsystem changes with respect to 
chemical potential at fixed @. We first focus on the second order transition case, which 
is presented in Figure (5). We can see that, from the insulator phase to the supercon- 
ductor phase, the entanglement entropy increases monotonically with the increment of 
chemical potential both in “deconfining phase” and “confining phase”. ° Furthermore, the 
entanglement entropy is continuous at critical chemical potential ue, but its slope has a 
discontinuous change at ue. Comparing with Figure (5), we can see a dramatic change 
in the first order case drawn in Figure (6). Although the entropy increases, there is an 
sudden jump in entropy as well as its slope at ue. The discontinuity or jump at critical 
point indicates some kind of significant reorganization of the degrees of freedom of the 
system, since new degrees of freedom are expected to emerge in the new phase. 

No matter of the order of the transition, the entanglement entropy increases as one 
increases chemical potential for the p-wave superconductor/insulator case. Such behavior 
is quite different from the s-wave superconductor/insulator case, where the entropy with 
respect to chemical potential is non-monotonic in the superconducting phase. However, 


Strictly speaking, the term phase transition here is inappropriate since the system itself, i.e., the state 
of the boundary field theory, does not change at all as one changes @. However, this behavior is reminiscent 
of many thermodynamic phase transitions in holographic calculations (see, for example Refs. [31, 32, 33]) 
and so here we adopt the terminology “phase transition” to convey this picture as in the literature. 

Due to the lake of numerical control at large chemical potential, we are not able to give the value of 
entropy for very large u. But up to u œ 15, we have confirmed that the trend of the entanglement entropy 
curve does not change. 
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Figure 6: The entanglement entropy as a function of chemical potential for strip with 
finite length along x direction. We choose a = 0.6 at fixed £ = 0.5L (within “deconfining 
phase”, left plot) and £ —+ oo (within “confining phase”, right plot), respectively. The 
transition from the insulator to superconductor is first order here. The critical chemical 
potential 4e ~ 2.157 is indicated by the vertical dotted blue line. In both plots, the solid 
curves come from the superconducting phase, while the dashed black lines come from the 
insulator phase. The physical curve is determined by choosing the dashed line below pe 
and solid upper curve above pe. 


we can find that the critical width as a function of chemical potential in the supercon- 
ducting phase behaves non-monotonic, which is presented in Figure (7). This behavior is 
qualitatively same as the s-wave case. 


3.2 Strip with finite width along y direction 


We now consider a straight belt geometry with a finite width £ along the y direction. 
The subsystem B locates in the slice r = + with the UV cutoff « — 0, and extends in x 
and 7 directions. The holographic dual surface yg extending to the bulk, is defined as a 
three-dimensional surface 


R R 
t=O, r =r) -7 <2 < 7 (R> o0), 0<n<T, (16) 


where R is the regularized length in the x direction. In this case, there also exist two kinds 
of dual surafces: connected configuration and disconnected configuration as in the previous 
subsection. Let us first concentrate on the connected configuration which is smooth at the 
turning point at r = r,. The dual surface yg starting from y = £ at r = I extends into the 
bulk until it reaches the turning point r = r,, then goes back to the AdS boundary r = 2 


at y = —£, 
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Figure 7: The critical width as a function of chemical potential in superconducting phase for 
a = 0.2 with ue œ 2.265 (left plot) and a = 0.6 with ue ~ 2.157 (right plot), respectively. 
For the first order phase transition case (right plot), the physical curve denoted by solid 
curve starts from the critical chemical potential pe ~ 2.157. 


In this case, the corresponding entanglement entropy is given by 


An L z 4 h(r)e-X) 2L 1 
oe = = T AP r g(r) (r)e eee A ee ee 2 ik RI(> te geny; 
K € 


(17) 
where the UV divergent part has been separated form the total entropy and S®” is a finite 
part. The width Z of the subsystem B and r, are related by the equation 


= -fah sf m (18) 
2 oe 2 (r) _r8g(r)h(r)e-x0) _ 1) 


Me Tgr) hraje x0 


Jex — regeare 


On the other hand, the disconnected configuration consists of two separated surfaces are 
located at y = +£, respectively, and extend into the bulk until they reach the tip r = ro, 
The entropy for this disconnected configuration is £ independent and given by 


isconnec 4nL z — xt) 27 L iscon 
Så = Sirr | rVJ/h(r)e~ 2 dr = ARG peer), (19) 
ro 


where once again, S4*’°°" denotes the finite part. 

We find that for the strip along the y direction, the entanglement entropy with respect 
to strip width @ behaves similar as the case along the x direction shown in Figure (4) 
for different choice of parameters, i.e., œ and u. As one tunes strip width, the “confine- 
ment /deconfinement” phase transition always exists. We present the entanglement entropy 
versus chemical potential in Figure (8) for a = 0.2 and in Figure (9) for a = 0.6, respec- 
tively. No matter of the order of the transition, the entanglement entropy increases as 
one increases chemical potential for the p-wave superconductor/insulator case. The crit- 
ical width as a function of chemical potential in the superconducting phase also behaves 
non-monotonic, which is presented in Figure (10). 
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Figure 8: The entanglement entropy as a function of chemical potential for strip with finite 
length along the y direction. We consider a = 0.2 at fixed @ = 0.5L (within “deconfining 
phase”, left plot) and £ —+ oo (within “confining phase”, right plot), respectively. The 
transition from the insulator to superconductor is second order here with the critical chem- 
ical potential ue ~ 2.265. In both plots, the solid curves come from the superconducting 
phase, while the dashed black lines come from the insulator phase. The physical curve is 
determined by choosing the dashed line below ue and solid upper curve above ue. 


Comparing the two setups, i.e., the strip with width Z along the x direction and the 
strip with width @ along the y direction, the results are qualitatively similar. However, 
the values of entanglement entropy and critical width in the former case are larger for 
given parameters than in the latter case and the difference grows quickly as one increases 
the chemical potential. Since the degree of anisotropy of the superconducting phase is 
characterized by the value of condensate. The growth of the difference of the entanglement 
entropy in two setups suggests that the entanglement entropy indeed describes the new 
degrees of freedom emerging in the superconducting phase. As a non-local quantity, the 
entanglement entropy for different configurations considered in this paper is not expected 
to have qualitatively difference, which is indeed confirmed by our numerical calculation. 
On the other hand, if considering local quantities for holographic p-wave model, such as 
conductivity [12, 16], one can find that the conductivity oyy orthogonal to the super current 
(ÎE) displays gapped dependence similar to the findings of s-wave case. The conductivity 
Ore parallel to the super current (J?) is qualitatively unlike oyy and the low-frequency 
behavior of Re[o,,| can be characterized very accurately in terms of the Drude model. The 
essential difference between two conductivities along different directions is the consequence 
of the anisotropy in the superconducting phase. 


4 Conclusion and discussions 


In our previous papers [28, 29], we studied the properties of entanglement entropy in the 
holographic s-wave superconductor /insulator transitions, where a non-monotonic behavior 
of the entanglement entropy with respect to chemical potential was found. Such behav- 
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Figure 9: The entanglement entropy as a function of chemical potential for the strip with 
finite length along the y direction. We consider a = 0.6 at fixed £ = 0.5L (within “decon- 
fining phase”, left plot) and l —>+ oo (within “confining phase”, right plot), respectively. 
The transition from the insulator to superconductor is first order here. The critical chem- 
ical potential ue ~ 2.157 is indicated by the vertical dotted blue line. In both plots, the 
solid curves come from the superconducting phase, while the dashed black lines come from 
the insulator phase. The physical curve is determined by choosing the dashed line below 
Hc and solid upper curve above pe. 


ior exists for more complex interaction [29]. In order to see whether the non-monotonic 
behavior is universal or not, in this paper we generalized our study to the p-wave case. 
Depending on the strength of the back reaction, the transition between the insulator phase 
and superconducting phase can be either second order or first order. However, unlike the 
case in s-wave model, there is no additional first order transition inside the superconducting 
phase. 

With the aim of providing more information on some universal properties of holographic 
superconductors which may shed toy-model-based insights on understanding the real high 
T, superconductors, we continued the study of the behavior of entanglement entropy in the 
holographic superconducting phase transition. We limited ourselves in this paper to inves- 
tigate an Einstein- Yang-Mills model, which is dual to the p-wave superconductor /insulator 
phase transition. In this model, the rotational symmetry acting on x and y is also broken 
in the superconducting phase. In order to account for such anisotropy, we calculated the 
entanglement entropy for two kinds of strip geometry with finite width @ along x direction 
and y direction, respectively. In both cases, the “confinement /deconfinement” phase tran- 
sition always appears. No matter of the order of the phase transition, the entanglement 
entropy as a function of chemical potential monotonically increases, which is much differ- 
ent from the s-wave case. Since the entanglement entropy is a non-local quantity in the 
field theory side and measures the area of a minimal surface extending into the bulk in 
the holographic description, we do not expect that the entanglement entropy for the two 
kinds of strip geometry will be qualitatively different. Our calculations indeed confirmed 
this and the entanglement entropy has a difference only quantitatively. More precisely, the 
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Figure 10: The critical width as a function of chemical potential in superconducting phase 
for a = 0.2 with ue ~ 2.265 (left plot) and a = 0.6 with ue œ 2.157 (right plot), respectively. 
For the first order phase transition case (right plot), the physical curve denoted by solid 
curve starts from the critical chemical potential ue ~ 2.157. 


difference increases as the condensate increases. Thus, it is reasonable to conclude that 
the entanglement entropy indeed describes the new degrees of freedom emerging in the su- 
perconducting phase. We also extracted the critical width which behaves non-monotonic 
with respect to chemical potential. This behavior is the same as the one in the s-wave 
case. In addition, we further provided the evidence that the entanglement entropy can in- 
deed indicate the appearance and the order of phase transition: In the second order phase 
transition case, the entanglement entropy is continuous, but its slope has a jump at the 
critical chemical potential, while in the first order transition case, both the entanglement 
entropy and its slope have a jump at the critical chemical potential. 

The different behavior of entanglement entropy in the s-wave case and p-wave case 
makes the aim to understand the non-monotonic behavior difficult at the moment. Nev- 
ertheless, after a second thought, the different behavior of entanglement entropy in both 
cases looks not so surprised, since the properties of holographic s-wave model and p- 
wave model are quite different. Ref. [22] studied (3 + 1)-dimensional s-wave superconduc- 
tor/conductor/insulator phase transition. For the conductor/superconductor case, the or- 
der of the phase transition is always second order. In contrast to the conductor /superconductor 
case, the insulator /superconductor phase transition is second order for small back reaction, 
for intermediate back reaction, there is an additional first order transition inside the su- 
perconducting phase and the transition finally becomes first order for sufficiently large 
back reaction. Similar phenomena exist even when one introduces more complex interac- 
tion [29]. Actually, we have also checked the same model in one dimension less, i.e., in 
(2 + 1)-dimensional case and found the same result. Let us now consider holographic p- 
wave models. The full back reacted p-wave superconductor/conductor transition in (3+ 1) 
dimensions studied in Ref. [15] shows that the order of the phase transition changes from 
second order to first order as one increases the strength of the back reaction, while for the 
(2+1)-dimensional case, the transition is always second order [27]. The (3+ 1)-dimensional 
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p-wave superconductor/insulator model studied in this paper and Ref. [23] behaves very 
similar to the (3 + 1) dimensional p-wave superconductor/conductor model. Let us note 
that the order of the phase transition for the p-wave case depends on not only the dimen- 
sion of space-time, but also the strength of the back reaction. This is quite different from 
the s-wave case. Furthermore, let us notice that in Ref. [1], the entanglement entropy for a 
(2 + 1)-dimensional chiral p-wave superconductor model is calculated by dividing the sys- 
tem as two parts separated by a hyperplane. In this model, there are four phases separated 
by three quantum critical points at chemical potential 4 = 0 and +4, which are labeled 
by the Chern number Ch as Ch = 0 (|u| > 4), Ch = —1 (—4 < u < 0), and Ch = +1 
(0 < u < +4). The entanglement entropy as a function of chemical potential is presented 
in Figure 4-(e) of Ref. [1]. One can see from the plot that the behavior of the entanglement 
entropy depends on not only the phase, but also the aspect ratio r = N,/N,, the ratio of 
lattice number in y direction over in x direction. In particular, in some phases both the 
monotonic and non-monotonic behaviors appear depending on the value of r. For a large 
r, the entanglement entropy is monotonic in a fixed phase while it becomes non-monotonic 
for a small r. It indicates that the entanglement entropy can indeed reveal many properties 
of many-body systems. 

In our setup, the holographic superconductor extends infinitely along x and y directions. 
In some sense this case corresponds to the large r situation discussed in Ref.[{1]. Therefore 
consider the entanglement entropy as a measure of the degrees of freedom for a given 
system, the monotonic behavior of the entanglement entropy in the holographic p-wave 
superconductor /insulator arises naturally. The normal insulator phase can be thought of as 
having little degrees of freedom due to the mass gap, which is consistent with the vanishing 
charge density. As we increase the chemical potential, the superconducting phase appears. 
The growth of the condensate and the charge density indicates that the total degrees 
of freedom increase, which just accounts for the increase of the entanglement entropy. 
However, it is still obviously necessary to understand the different behavior of entanglement 
entropy in holographic s-wave and p-wave superconductor/insulator transitions. For this, 
one needs know the details of dual field theories for these two kinds of models. In principle, 
one can finish this task if one could embed these two models into string theory or M- 
theory. Unfortunately, we have not yet had the microscopic theories for these two models. 
Nevertheless, the phenomenological bottom-up approach we adopted here can provide some 
basic information of a particular class of strongly coupled field theory. In this sense it might 
also helpful to use the entanglement entropy as a probe to study the holographic d-wave 
models [34]. 
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A Numerical details 


In this Appendix, we give the details for solving the coupled equations of motion (7). In 
order to find the solutions for all the six functions F = {¢,w, f, h, g, X}, one must impose 
suitable boundary conditions at boundary r — oo and at the tip r = rọ. The general falloff 
near the boundary is given by (8) where we choose wo = 0 since one wants the condensate 
to arise spontaneously. 

Regularity of the solution at the tip r = ro requires that all six functions have finite 
values at the tip and have Taylor series expansions near the tip 


F = Fire) + Mo) — 10) ees (20) 


Plugging the expansion (20) into (7), and using g(ro) = 0, we are left with six independent 
parameters {79, @(ro), (ro), f (ro), A(ro), x(ro)}. However, the equations of motion (7) have 
four useful scaling symmetries 


f>Xf, oS Ne, (21) 

hN h whe, (22) 

L>x+à, 9 ne? (23) 

rar, (t,2,y,n) 2A *(E2,y,n), (¢,w) > A(b,w). (24) 


Taking advantage of such four scaling symmetries, we can firstly set {rọ = 1, f(ro) = 
1,h(ro0) = 1,x(r0) = 0} for performing numerics. After solving the coupled differential 
equations, we should use the first three symmetries again to satisfy the asymptotic condi- 
tions. We will choose $(r9) as a shooting parameter to match the source free condition, 
i.e, Wo = 0. 

All in all, for each choice of œ and w(ro), we can solve the equations of motion (7). 
After solving the equations, we can obtain the current condensate (J?) = 20" we, chemical 
potential u = o and charge density p = 20° be by reading off the coefficients w2, do and 
@2 in (8), respectively. 

In fact, it is convenient in the numerical calculations to make a coordinate transfor- 
mation from r coordinate to z coordinate by defining z = 1/r. In that case, the infinite 
boundary is now at z = 0 and the tip sits at z9 = 1/rp = 1. Two typical sets of solution 
for the metric and gauge field configurations are presented in Figure (11) for œ = 0.2 and 
a = 0.6. respectively. Note that the metric functions g and h will be used in calculating 
the holographic entanglement entropy. For fixed strength of back reaction a, one has a one 
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Figure 11: Typical soliton solution with nonvanishing vector hair. The solid curves stand 
for the case with a = 0.6, where the value of w at tip is w(zọ) ~ 1.576 and the corresponding 
identification length I œ 2.842, while the dashed lines are for the case with a = 0.2, where 
the value of w at tip w(zo) ~ 1.922 and the corresponding identification length F ~ 3.118. 


parameter family of solutions denoted by the value of w at the tip. However, for different 
choices of w(ro), the identification length I in 7 coordinate will be different. In order 
to compare different solutions appropriately, the boundary geometry should be the same. 
Making use of the scaling symmetry (24), one can find that the relevant quantities scale 
as follows 


= T/A, ro Aro, H =>> AL, p> *p, ip > AT 
gcon/ discon = BS ae aaa l> lJA, 


Scon/discon 


(25) 


where stands for the finite part of the entanglement entropy defined in section (3) 
and £ is the strip width. Since in our units setup, the identification length I in the pure 
soliton is 7, we will scale all for each solution to be m. It should be stressed here that 
after the scaling transformation, the tip rọ will be no longer at rọ = 1. 

For each a, as one increases the chemical potential, the solution with non-vanishing 
vector “hair” will appear and will be finally thermodynamically favored over the pure AdS 
soliton solution above the critical chemical potential ue, where a phase transition from 
insulator phase to superconducting phase happens. We present critical chemical potential 
Uc as a function of the strength of the back reaction a in Figure (12). For small a, the 
transition is second order and ju, is almost independent of the strength of the back reaction. 
For larger a, the transition becomes first order and the critical chemical potential decreases 
as a increases. Therefore, there exists a critical value a, ~ 0.538 + 0.002, beyond which 
the phase transition changes from second order to first order. 
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Figure 12: The critical chemical potential u« as a function of the strength of back reaction 
a. The critical strength a, ~ 0.538 + 0.002, denoted by the dashed vertical line, where the 
order of transition is changed. Here we scale all I’ for each a to be 7. 
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